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A covariant gauge independent derivation of the generalized dispersion relation of electromag- 
netic waves in a medium with local and linear constitutive law is presented. A generalized photon 
propagator is derived. For Maxwell constitutive tensor, the standard light cone structure and the 
standard Feynman propagator are reinstated. 



i> : 
o : 
o ■ 

(N ■ 



I 

Oh' 



> 

p 

o 
i> 
o 



X 



PACS numbers: 04.20.Cv, 04.50.+h, 03.50.De 

I. INTRODUCTION 

From a formal point of view [1] , [2] , the Maxwell elec- 
trodynamics theory can be represented by a system of 
two independent equations 



= J' 



(1.1) 



where two independent antisymmetric tensors, the field 
strength tensor Fij and the excitation tensor density i/*-' 
are involved. The electromagnetic current vector density 
is denoted by J* Here the commas stand for ordinary 
derivatives, the indices range from to 3, the Levi-Civitas 
permutation tensor is normalized by e^^^^ _ ^ 

For most applications it is enough to assume a local, 
linear, homogeneous constitutive relation between the 



fields F,j and H'^ , 



kl ■ 



(1.2) 



By the definition, the constitutive tensor has to 

respect the symmetries of the fields Fij and H^^ , 



X 



ijkl 



^[ij]kl _ ^ij[kl] 



(1.3) 



Hence it has, in general, 36 independent components. 

The standard Maxwell electrodynamics in vacuum is 
reinstated in this formalism by a special choice of the 
Maxwell-Lorentz constitutive tensor 



(Max)^ijfci 



g^'g'") ■ 



(1.4) 



Here g^^ is the Lorentz metric components, while Aq is a 
constant with the dimension of an admittance. 

In this paper, we study the formal scheme (1.1), (1-2) 
with a general constitutive tensor The physical 

spacetime is considered as a bare manifold without met- 
rics or connection. All the information about the geom- 
etry of this space is encoded in the constitutive tensor. 
In other words we are dealing with a premetric electro- 
dynamics. 

Such a construction is applicable for description of a 
rather wide range of physics effects. As a classical field 
theory, the premetric electrodynamics involves the stan- 
dard Maxwell electrodynamics in vacuum and even pro- 
vide a possibility to describe the additional degrees of 



freedom (axion, dilaton and skewon) as the premetric 
partners of photon [3], [4]. Moreover, since the metric is 
a secondary quantity in this scheme, its form [5], [6], [7], 
[8] , and the signature [9] are derived from the properties 
of the constitutive tensor. The nonminimal coupling of 
the electromagnetic field to the torsion yields the bire- 
fringence of vacuum [10], [11]. This effect find its natural 
description in the premetric scheme, [12], [13]. 

Another interesting area of application is the models 
with violation of Lorenz invariance. In particular the 
Carroll-Field- Jackiw modification of the Maxwell electro- 
dynamics [14], see also [15], is embedded in the premetric 
scheme. The wave propagation in this model requires, 
however, to go beyond the geometrical optics approxi- 
mation [16]. This problem will be considered in a con- 
tributed publication. 

The mathematical methods similar to used here was 
shown to be useful in ray optics applications to GR [17] 
and in quantum plasmadynamics [18]. 

In the present letter, we give a covariant gauge inde- 
pendent derivation of the generalized dispersion relation 
for the premetric electrodynamics. Moreover, we derive 
a generalized Green function in the momentum represen- 
tation ~ a generalized photon propagator. 



II. DISPERSION RELATION 

To study the wave propagation in the premetric elec- 
trodynamics model, we solve the first equation of (1.1) in 
term of potentials Fij — {l/2){Aij — Aj^i). Substituting 
it into (1.2) and (1.1) and the current J* to be equal to 
zero we derive 



X 



ijkl 



k.l] 



0. 



(2.1) 



To study the wave-type solutions of this equation we con- 
sider an anzatz 



Aij [x) 



(2.2) 



where tp = ip{x'') while is a constant covector. Such 
solutions always exist on sufficiently small neighborhoods 
even on bare manifold [17]. Denote the wave covector as 
9i = V.i- 



2 



In the geometrical optics approximation, the changes 
of the media parameters are neglected relative to the 
changes of the wave characteristics. Consequently we 
come to an algebraic system 



M'''ak = , where M''' = ]-x'^^^ 



X' QiQj 



(2.3) 



Due to the symmetries of the constitutive tensor (1.3) 
the matrix of the system satisfies 



M^'^Qk = , M^'^Qi = . 



(2.4) 



The first relation of (2.4) means the gauge freedom of the 
vector potential while the second relation is interpreted 
as a charge conservation condition. Due to (2.4), the 
rows (and the columns) of the matrix M'-' are linearly 
dependent, so its determinant is equal to zero. Moreover, 
the gauge relation (2.4) can be interpreted as a fact that 



Ok = Cqk 



(2.5) 



is a formal solution of (2.3). This solution does not give 
a contribution to the electromagnetic field strength so it 
is unphysical. 

An additional physically meaningful solution has to be 
linear independent on (2.5). A linear system has two or 
more linear independent solutions if and only if its rank 
is two (or less). Consequently, a generalized electrody- 
namics system has a physically meaningful solution if 



Aij = . 



(2.6) 



Here we involved the adjoint matrix Aij - a matrix con- 
structed from the cofactors of M'^ . The components of 
the adjoint matrix are expressed by the derivatives of the 
determinant relative to the entries of the matrix 



_ ddetjM) 



1 
3! 



-Jlll2»3"^JJlj2j3 



(2.7) 

Since the adjoint matrix has, in general, 16 independent 
components it seems that we have to require 16 indepen- 
dent conditions. The following algebraic fact shows that 
the situation is rather simpler. 

Proposition: If a square n x n matrix M*-' satisfies 
the relations 



(2.8) 



for some nonzero vector qi, its adjoint matrix Aij is rep- 
resented by 



Aij = \{q)qiqj ■ 



(2.9) 



For a formal proof of this fact, see [20]. Consequently, 
instead of (2.6), we have only one condition 



X{q)=0. 



(2.10) 



This condition is necessary to have physically meaning- 
ful solutions of the generalized wave equation, so it is a 
generalized dispersion relation. 



The problem now is to derive from (2.9) the explicit 
expression for the function X{q). It is provided [20] by 
using the fact that the functions involved in (2.9) are 
homogeneous polynomials. In fact, A^j is of the sixth 
order in the wave covector g*, while A(<3') is of the fourth 
order. Applying twice the derivatives with respect to the 
components of the covector and using Euler's rule for 
the homogeneous functions, we obtain 



A(g) = 



1 d^Ai 



72 dqidqj 



(2.11) 



In term of the matrix M*-' , the function X(q) is rewrit- 
ten as 



dqi dqj 



^M'^^K (2.12) 



This expression may be useful for actual calculations of 
the dispersion relation for different media. 

In order to have an explicit expression of the function 
A in term of the constitutive tensor we have to calculate 
the corresponding derivatives. The resulting dispersion 
relation is 

p.. . . f.. . . ( ^n(v)3i ^i2ab]2 , 

tjlil2l3tjjlj2j3 I A A I 

4x^^^*"^^'^x*^^^'^^'^)x^^^*^«ag6«cgd = 0. (2.13) 

This equation is completely equivalent to the recently 
proposed [2] covariant dispersion relation 

e«ii2i3eiiu2i3X"^^"x''^^'^'X*^''=^'=*9a965c% = . (2.14) 

Indeed, in the special coordinate basis with qi = 
(q, 0,0, 0), both equations yield the same non-covariant 
expression. Also the direct proof of the equivalence of 
two forms is acceptable [19]. 

The function A(g) is a fourth order polynomial. When 
it is separated to a product of two non-positive defined 
quadratic factors birefringence effect emerges. This effect 
is well known from the classical optics. However, in the 
premetric approach two light cones explicitly represent 
violation of Lorentz invariance. The non-birefringence 
condition can be given [7] in a rather simple covariant 
form: For an arbitrary covector q, 



X{q) > 



(2.15) 



has to be satisfied. For a component-wise representation 
of this condition, see also [6]. 

Observe an important special casc\ When the skewon 
part absents, the constitutive tensor respects the symme- 
tries 



(2.16) 
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In this case, two terms in (2.13) are proportional one 
to another. Thus two additional expressions for the re- 
stricted dispersion relation emerge 

tiiii2»3'^JJlj2j3 A A A 1a1b<dc<id — 'J , K'^--^') 

and 

Cjtit2i3"^JJij2j3 A A A Halb^c'id — 

(2.18) 

For the Maxwell constitutive tensor (1.4), the matrix 

M*-' takes the form 

("^=^)M'^ = AoT^ - q\^) . (2.19) 

The corresponding adjoint matrix is 

= - {\osF^\\qj . (2.20) 

Consequently, in this special case, the dispersion relation 
takes its regular form =Q. 



III. PHOTON PROPAGATOR 

Let us return to the full inhomogeneous Maxwell equa- 
tion with a non-zero current. In the "momentum" repre- 
sentation, it takes the form 



(3.1) 



Observe that the charge conservation law is expressed 
now as 



(3.2) 



It is useful to have a formal solution of the equation (3.1) 
for a an arbitrary given current jf. . Such a solution is usu- 
ally given by the Green function or photon propagator, 
Dij{q). This tensor is defined in such a way that the 
covector 



ak = -Dkif 



(3.3) 



is a formal solution of (3.1). Note that, due to the 
gauge invariant and charge conservation, the propagator, 
Dij{q), is defined only up to addition of terms propor- 
tional to the wave covector qi, 



Dij Dij + (j)iqj + ipjqi . 



(3.4) 



Here the components of the covectors 0i and tpi are 
arbitrary functions of the wave covector. In the stan- 
dard electrodynamics, an expressions for this quantity is 
known as the Feynman propagator 



9^J 



(3.5) 



Note that it is expressed by a symmetric matrix. Thus 
also the covectors (pi and ipi are usually taken to be equal 



one to another. In our general setting, Dij can be asym- 
metric. Consequently it useful to preserve two arbitrary 
covectors in (3.4). 

Substituting (3.3) into (3.1) we get 



{M'>'Dkm-Sl)r = 0, 



(3.6) 



Note that the matrix M'*^ is singular, so the propagator 
cannot be taken to be proportional to the inverse of M*'' . 
In the standard case, it is not a problem. Indeed, with 
the matrix M' given by (2.19) and with the Feynman 
propagator (3.5) we have 



Aq 



q' 



-9- 



(3.7) 



When this expression is substituted in (3.6) the second 
term disappears due to the charge conservation equation 
(3.2) and the equation is satisfied. 

Our task is to derive an expression for the photon prop- 
agator in a general case when the metric tensor is not 
acceptable. To deal with the singular matrix M'^, we 
consider the tensor density 



ijkl 



adet(M) 



(3.8) 



This is the, so called, second adjoint (or the second ad- 
jugate compound) of the matrix IvP^ . It is obtained by 
removing two arbitrary rows and two arbitrary columns 
from the original matrix. Observe that due to its defini- 
tion the second adjoint tensor respects the symmetry 

Bijki = Bkiij (3.9) 

It is expressed by the components of the matrix M*^ as 

(3.10) 



1 



Bijki = ^eiki,i,ejij,j,M^^^^M^^^^. 



From this expression, we read ofi^ the additional symme- 
tries 



B 



ijkl 



-B, 



kjil 



-B 



ilkj 



(3.11) 



Let us derive an identity involving the second adjoint 
tensor. The derivative of the generalized Laplace expan- 
sion AijM'^'' = relative to the enters of the matrix Af* 
yields 



BijrsM 



ik 



-A ■^'^ 



(3.12) 



We multiply now both sides of the equation (3.1) by 



the tensor B,jrs to get 
Using (3.12) we rewrite it as 



-Bi 



(3.13) 



(3.14) 
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Substituting (2.9) we get 



In term of the matrix M '-^ it takes the form 



(3.15) 



We are coming once more to the same problem: How to 
"divide" both sides of this equation by the covcctor qi in 
a covariant manner? Observe that A and Birnnk are ho- 
mogeneous polynomials in q of the order 4. Assuming 
to be independent on q, we see that is a homogeneous 
polynomial in q of the order —2. Note that this is in a 
correspondence with the classical expressions (3.5). 

Applying twice the partial derivatives with respect to 
the components of the wave covector and using Euler's 
rule for the homogeneous functions we come to 



1 ^2 



6 dqmdq„ 



B. 



imnk 



(3.16) 



Consequently we derived an expression for the gener- 
alized photon propagator 



52 



(3.19) 

And finally we derive an expression the generalized pho- 
ton propagator via the constitutive tensor 



84A 



^ji (mn)mi ^j2abTn2 



_|_2-^Jl('"«)"»l-^i2("*)'"2 



QaQb ■ 



(3.20) 



For the Maxwell constitutive tensor, the second adjoint 
takes the form 

Bijki = '^>'lgq^ {gijqm +9kmqj) - {guqjqk + gkjqiqi) ■ 

(3.21) 

Calculating with (3.18) we come to the standard Feyn- 
man propagator expression. 



1 92 



6( 



%dqn 



B 



rrnjn 



(3.17) 



Using the homogeneity of the polynomials involved here 
we get certain equivalent expressions 



Di 



1 



92 



1 d^Br, 



42A dqmdqn 42A dqmdqr 
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